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Abstract 

A narrow resonance in D* ~p and D* + p invariant mass combinations is observed in in- 
elastic electron-proton collisions at centre-of-mass energies of 300 GeV and 320 GeV at 
HERA. The resonance has a mass of 3099 ± 3 (stat.) ± 5 (syst.) MeV and a measured 
Gaussian width of 12 ± 3 (stat.) MeV, compatible with the experimental resolution. The 
resonance is interpreted as an anti-charmed baryon with a minimal constituent quark com- 
position of uuddc, together with the charge conjugate. 
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1 Introduction 



Several experiments have recently reported the observation of a narrow resonance with mass 
in the region of 1540 MeV, decaying to K + n or K®p [1]. This state has both baryon number 
and strangeness of +1, such that its minimal composition in the constituent quark model is 
uudds. It has thus been interpreted as a pentaquark [2, 3], the 9 + . There is also evidence for 
two related states with strangeness of —2 [4]. Various models have been put forward to explain 
the nature of these states and the structure of the multiplet that contains them [3, 5, 6]. The 
possibility of pentaquark states in the charm sector has also been discussed [7], with renewed 
theoretical interest in calculating their expected properties [5, 8] following the observation of 
strange pentaquarks. 

The electron-proton collider, HERA, is a copious producer of both charm and anti-charm 
quarks, the dominant production mechanism being boson-gluon fusion, — > cc. The spec- 
troscopy of charmed hadrons can be studied using the final states to which the quarks and 
anti-quarks hadronise. This paper reports the first evidence for a baryon with exotic quantum 
numbers in the charm sector, using deep inelastic scattering (DIS) data taken with the HI de- 
tector. A resonance is observed when combining D*~ — > D°n~ — > K + n~nJ candidates with 
proton candidates and when combining D* + — > -D°7T+ — > K~it + it+ candidates with antiproton 
candidates. 1 The resonance is also observed in an independent photoproduction data sample. 



2 Experimental Procedure 
2.1 HI Apparatus 

The tracks from charged particles used in this analysis are reconstructed in the HI central 
tracker, whose main components are two cylindrical drift chambers, the inner and outer cen- 
tral jet chambers (CJCs), covering the polar angle region 2 20° < 9 < 160°. The inner and 
outer CJCs are mounted concentrically around the beam-line, have 24 and 32 sense wires, re- 
spectively, and cover radii between 20 cm and 84 cm. The information from the CJC sense 
wires is digitised using 100 MHz FADCs, providing simultaneous charge and timing measure- 
ments. The CJCs lie within a homogeneous magnetic field of 1.15 T, which allows measure- 
ments of the transverse momenta of charged particles. Two additional drift chambers com- 
plement the CJCs by precisely measuring the z coordinates of track segments and hence as- 
sist in the determination of the particle's polar angle. The Central Silicon Tracker, consisting 
of two layers at radii of 6 cm and 10 cm, is also used to improve the charged particle track 
and event vertex reconstruction. The transverse momentum resolution of the central tracker is 
a(p T )/p T ~ 0.005 p T [GeV] ©0.015. The charge misidentification probability is negligible 

'in the remainder of this paper, particle charges are not generally given. Both charge conjugate configurations 
are always implied, unless explicitly stated otherwise. The notation n s is used to distinguish the low momentum 
pion released in the D* decay from that from the D° decay. 

2 The HI experiment uses a coordinate system in which the positive z-axis is defined by the direction of the out- 
going proton beam. The polar angle 9 of a particle is defined relative to this axis and is related to the pseudorapidity 
rj by r\ = — In tan 9/2. 
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for particles originating from the primary vertex which have transverse momenta in the range 
relevant to this analysis. 

The specific ionisation energy loss of charged particles is derived from the mean of the 
inverse square-root of the charge collected by all CJC sense wires with a signal above threshold. 
The resolution is a (dE/dx) /(dE/dx) ~ 8% on average for minimum ionising particles [9]. 

A lead/scintillating-fibre spaghetti calorimeter (SpaCal) is located in the direction of the 
outgoing electron beam. It contains both electromagnetic and hadronic sections and is used to 
detect the scattered electron in DIS events. The global properties of the hadronic final state are 
reconstructed using an algorithm which takes information from the central tracker, the SpaCal, 
and also from a Liquid Argon calorimeter, which surrounds the central tracker. The DIS events 
studied in this paper are triggered on the basis of a scattered electron in the SpaCal, comple- 
mented by the signals in the CJCs and multi-wire proportional chambers in the central tracker. 
Further details of the HI detector can be found in [10]. 

2.2 The DIS Data Sample 

The analysis is carried out using data taken in the years 1996-2000, when HERA collided elec- 
trons 3 of energy 27.6 GeV with protons at 820 GeV (1996-1997) and 920 GeV (1998-2000). 
The integrated luminosity of the sample is 75 pb _1 . 

The scattered electron energy, measured in the SpaCal, is required to be above 8 GeV, and 
the virtuality of the exchanged photon 4 is required to lie in the range 1 < Q 2 < 100 GeV 2 , as 
reconstructed from the energy and polar angle of the electron. To ensure that the hadronic final 
state lies in the central region of the detector, the inelasticity of the event is required to satisfy 
0.05 < y < 0.7, calculated using the scattered electron kinematics. The z coordinate of the 
event vertex, reconstructed using the central tracker, is required to lie within 35 cm of the mean 
position for ep interactions. The difference between the total energy E and the longitudinal 
component of the total momentum p z , calculated from the electron and the hadronic final state, 
is restricted to E — p z > 35 GeV. This requirement suppresses photoproduction background, 
where a hadron fakes the electron signature. 

2.3 Selection of D* Meson and Proton Candidates 

The decay channel D* — > D°n s — > Kitit s is used to reconstruct D* mesons. The charged 
particle selection criteria, dE/dx requirements and transverse momentum cuts on the decay 
products are very similar to those used in previous HI analyses [11]. Unlike-charge particle 
combinations are made to form K t tt ± pairs, where the particles are required to be consistent 
with kaons and pions according to their dE/dx measurements. Those combinations that give 
rise to an invariant mass within 60 MeV of the nominal D° mass of 1864.5 MeV [12] are 

3 The analysis uses data from periods when the beam lepton was either a positron (88% of the total) or an 
electron (12% of the total). 

4 The inclusive DIS kinematic variables are defined as Q 2 = —q 2 , y = q ■ p / k ■ p and x = —q 2 / 2q ■ p, where 
q, k and p are the 4-vectors of the exchanged photon, the incident electron and the incident proton, respectively. 
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combined with further pion candidates (n s ) with opposite charge to the kaon. To obtain good 
experimental resolution and background rejection in the D* reconstruction, the standard mass 
difference technique [13] is used, based on the variable 

AM D * = m(K7vn s ) - m(K7i) . (1) 

The sample is restricted to a region where the background to the D* signal is relatively small by 
requiring that the D* candidates have transverse momentum p T (D* ) > 1 .5 GeV, pseudorapidity 
— 1.5 < i](D*) < 1 and production elasticity z(D*) — (E — p z ) D */2yE e > 0.2, where E e is 
the electron beam energy. 

The resulting AM D . distribution is shown in figure la. Here and elsewhere in this paper, the 
error bars shown represent the square-root of the numbers of entries in each bin. A prominent 
signal on a smooth background is observed around the expected D* — D° mass difference. 
The AM D * distribution is also shown in figure la for a "wrong charge D" background sample, 
which is formed by K ± n ± combinations in the accepted D° mass range. The "wrong charge 
D" distribution gives a good description of the correct-charge D° combinations away from the 
D* peak. Candidate D* mesons for which AM D * lies within ±2.5 MeV of the nominal mass 
difference m(D*) — m(D°) = 145.4 MeV [12] are selected for further analysis. 

Proton candidates are selected using requirements on the particle dE/dx measurements. 
Figure lb shows the dE/dx values, plotted against momentum, for a sample of particles which 
yield a mass M(D*p) < 3.9 GeV (see equation 2) when combined with Kirir s candidates 
falling in the accepted AM D * region. The likelihoods that a particle is a pion, kaon or pro- 
ton are obtained from the proximity of the measured dE/dx to the most probable values for 
each particle type at the reconstructed momentum. The most probable dE/dx values are de- 
rived from phenomenological parameterisations [9], shown in figure lb, which are based on the 
Bethe-Bloch formula. The normalised proton likelihood L p is defined to be the ratio of the pro- 
ton likelihood to the sum of the pion, kaon and proton likelihoods. For momenta p(p) < 2 GeV, 
a requirement L p > 0.3 is applied, which selects protons where they are clearly identified at 
low momentum and suppresses contributions close to the crossing points of the proton, pion 
and kaon parameterisations. For p(p) > 2 GeV, the requirement is loosened to L p > 0.1, which 
suppresses background from particles with large dE/dx such as electrons. The main selection 
criteria are summarised in table 1 . 



3 Analysis of D*p Combinations 
3.1 D*p Invariant Mass Distributions 

The D* and proton candidates are combined to form the mass difference m(Kinr s p)—m(KiTiT s ), 
to which the D* mass of 2010.0 MeV [12] is added to obtain the mass of the D*p combination. 
The distributions in 

M(D*p) = m{Kmr s p) - m(Kmi s ) + m(D*) (2) 

are shown in figure 2a for "opposite-charge" D*p combinations (K~7v + irfp and K + n~n~p) 
and in figure 2b for "same-charge" D*p combinations (K~n + iT+p and K + it~it~p). A clear 
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Figure 1: (a) AMp* distribution for K T 7i ± 7rf combinations as described in the text. For com- 
parison, the distribution from "wrong charge D" combinations, where the K and n yielding the 
D° mass have the same charge, is also shown, (b) Specific ionisation energy loss relative to 
that of a minimally ionising particle, plotted against momentum, for the sample described in 
the text. The curves indicate parameterisations of the most probable responses of the CJCs for 
pions, kaons and protons. 
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D° 


p T (K) > 500 MeV 
p T (n) > 250 MeV 
p T (K) +p T (n) > 2 GeV 
\m{Kir) - m(D°)\ < 60 MeV 


D* 


p T (ir s ) > 120 MeV 

\AM D * - m(D*) + m(D°)\ < 2.5 MeV 
p T (D*) > 1.5 GeV 
-1.5 < rj(D*) < 1 
z(D*) > 0.2 


P 


p T (p) > 120 MeV 

L p > 0.3 for p(p) < 2 GeV 

L p > 0.1 foe p(p) > 2 GeV 



Table 1 : Summary of the kinematic and proton energy loss selection criteria applied to define 
the D* and proton candidates. 



and narrow peak is observed for the opposite-charge combinations at M(D*p) ~ 3100 MeV. 
Approximately half of the events in this signal arise from each of the D*~p and D* + p com- 
binations (see section 4). The distribution for the same-charge combinations shows a small 
enhancement in the M(D*p) region in which the opposite-charge signal is observed. 

The background distributions for the D*p combinations are modelled by the sum of two 
contributions, which are shown in figure 2. Background from random combinations not involv- 
ing charm is modelled using the "wrong charge D" combinations, as described in section 2.3, 
combined with proton candidates as for the correct-charge D° sample. Combinatorial back- 
ground from D* mesons with real or misidentified protons is modelled using simulated events 
from the RAPGAP [14] Monte Carlo model applied to D* production in DIS, including string 
fragmentation and decays from JETSET [15, 16]. The RAPGAP model gives a good descrip- 
tion of the shapes of the inclusive D* distributions. This contribution ("D* MC" in figure 2) is 
normalised according to the D* yield in the data (figure la). 

No significant structures are observed in either component of this background model. The 
model gives a reasonable description of the shape and normalisation of the M(D*p) distribution 
away from the signal region for the opposite-charge combinations. The M(D*p) distribution 
from the same-charge combinations is also well described in shape, though the model prediction 
lies approximately 15% above the data. 

Alternative models have been studied for the background distribution for the opposite- 
charge D*p combinations. Similar distributions to those shown in figure 2 are obtained when 
a DJANGO [17] simulation of inclusive DIS is used to replace both model components. The 
same is true when the RAPGAP model of the _D*-related background is replaced by simulations 
with modified parton shower dynamics (CASCADE [18]) or fragmentation (HERWIG [19]). 
In all cases, no resonant structures are observed in the simulated M(D*p) distributions. Pos- 
sible contributions from beauty decays have been considered using a further RAPGAP Monte 
Carlo simulation. After normalising to the luminosity of the data, the resulting contribution is 
negligible. 



8 




3 3.2 3.4 3.6 

M(D*p) [ GeV ] 




Figure 2: Distributions in M(D*p) for (a) opposite-charge and (b) same-charge D*p combina- 
tions. The data are compared with a two-component background model in which "wrong charge 
D" K ± n ± combinations are used to describe non-charm related background and the "D* MC" 
simulation describes background involving real D* mesons. 
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The events giving contributions in the peak region of the opposite-charge M(D*p) distribu- 
tion have been visually scanned and no anomalies are observed in the events or the candidate 
tracks. All entries within ±24 MeV of the peak arise from different events, with one exception 
where the same n, ir s and proton candidates contribute with two different K candidates. For 
the full M(D*p) range shown in figure 2a, there are an average of 1.12 entries per event. The 
signal is present in each case when the data are divided into two sub- samples of similar size, 
discriminated in variables such as x or Q 2 , the pseudorapidity or transverse momentum of the 
D*p composite, or the data taking period. The peak also remains clearly visible for reasonable 
variations in the binning or selection criteria. In all cases, the observed mass and width of the 
peak are stable to within a few MeV. 

3.2 Particle Identification Tests 

The D* and proton content of the signal in the mass distribution from opposite-charge D*p 
combinations has been investigated in complementary studies. The D* content is tested by 
forming the M(D*p) distribution (equation 2) with the full proton selection, but with no re- 
quirement on AM D *. The AMp. distribution (equation 1) is shown in figure 3 for events in a 
±15 MeV window around the D*p signal (3085 < M(D*p) < 3115 MeV). For comparison, 
a similar distribution is shown, taken from side bands with 2990 < M(D*p) < 3070 MeV and 
3130 < M{D*p) < 3210 MeV, scaled by a factor of 3/16 to account for the different widths 
of the sample regions. Away from the D* peak, the distribution in AM D . from the M(D*p) 
side bands gives a good description of that from the M(D*p) signal region, in both shape and 
normalisation. However, there is a clear difference around the expected value of AM^. for 
true D* mesons, where the distribution from the signal region lies well above that from the side 
bands. The signal region in M(D*p) is thus significantly richer in D* mesons than is the case 
elsewhere in the distribution. 

The protons are clearly identified at low momentum, where the most probable dE/dx for 
protons is well separated from those for other particle species. The analysis has been repeated 
with the proton momentum restricted to p(p) < 1.2 GeV, the measured dE/dx required to 
be larger than that for a minimum ionising particle by a factor of at least 1.15 and the proton 
likelihood requirement modified to L p > 0.5. With this tighter proton selection, the require- 
ments on z(D*) and p T (K) + p T (n) are removed. In figure 4, the M(D*p) distribution for this 
selection is compared with the predictions of the background model described in section 3.1. 
The enhancement in the region M(D*p) ~ 3100 MeV remains visible. The candidate proton 
tracks in the signal region (3085 < M(D*p) < 3115 MeV) have an average (L p ) = 0.92. 

Several further checks of the particle identification have been carried out using the data. No 
signal is obtained when the D* selection is modified such that the "wrong charge D" candidates 
are taken (figure 2a) or when the D* candidates are replaced by Knn s combinations from a side 
band of the AM D . distribution. No clear peaks are observed when the proton candidates are 
assumed to have the kaon or pion mass. The signal is also not observed when the proton mass 
hypothesis is retained, but the dE/dx requirement is modified to select pion or kaon candidates 
rather than proton candidates. 
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Figure 3: AMp* distribution for events in a 30 MeV window about the signal in the opposite- 
charge M(D*p) distribution, with no requirement on AM D «, compared with the corresponding 
distribution from side bands in the M(D*p) distribution, normalised according to the widths of 
the chosen sample regions. 



3.3 Kinematic and Reflection Tests 



Possible kinematic reflections that could fake the signal have been ruled out by studying in- 
variant mass distributions and correlations involving the K, n, ir s and proton candidates under 
various particle mass hypotheses. For example, there is no evidence for any resonant structure 
or correlations with M(D*p) in the invariant mass combinations m(Kp), mijip) or m(n s p) of 
the proton candidates with the decay products of the D* meson. 

Detailed studies have been carried out of the contribution to the M(D*p) distribution from 
the neutral, orbitally excited, P-wave D°(2420) and D\ *(2460) mesons and their charge conju- 
gates [12], both of which decay to D* ± n T . A simulation of the —> D*n and D®* — > D*n 
decays with the PYTHIA [20] Monte Carlo generator is used to estimate their contribution to 
the observed signal. The simulated widths are set to the results from recent measurements [21] 
and the normalisation is obtained from the observed D\ and D®* yields in the data, as obtained 
from the distribution in M(D*n) = m(K T 7r ± Tcf7r T )—m(K T 7r ± Tcf)+m(D*). The reflections 
due to the D\ and D® * mesons when the decay pion is misidentified as a proton yield a broad 
distribution in M(D*p), with a maximum below the signal region. The predicted contribution 
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Figure 4: M(D*p) distribution for a low momentum proton selection with p(p) < 1.2 GeV. 
The data are compared with a two-component background model in which "wrong charge D" 
K ± ir ± combinations are used to describe non-charm related background and the ll D* MC" 
simulation describes background involving real D* mesons. 



within ±24 MeV of the observed peak is approximately four events. The signal in the measured 
M(D*p) distribution covers the full available phase space in M(D*ti), with no evidence for en- 
hancements in the region of the D\ and D®* mesons. The possibility of reflections involving 
the D s i(2536) or _D s j (2573), are similarly ruled out. 5 

The kinematics of the D* and proton candidates from the decay of a resonance would be 
expected to be different from those of the background distribution. Such a difference is observed 
for the opposite-charge D*p signal, as illustrated in figure 5a. The momentum distribution p(p) 
is shown for all particles of opposite charge to the D* candidate which lead to entries in the 
signal and side band regions of M(D*p), as defined in section 3.2. No requirements are placed 
on the proton likelihood. The two side bands with larger and smaller M(D*p) than the signal 
give rise to compatible momentum spectra, both of which are significantly softer than that in 
the signal region. A similar difference is observed when the M(D*p) side band is replaced by 
combinations which lie in the signal region of M(D*p), but fall in a side band of the AMp* 
distribution. 

Figure 5a suggests that, with no proton dE/dx requirements, the signal-to-background ratio 
improves as p(p) increases. In figure 5b, the M(D*p) distribution is shown for p(p) > 2 GeV, 

5 Given the proximity of the mass of the observed resonance to the J/ ip mass, possible backgrounds involving 
J I ip decays have been considered. Baryon number and other conservation laws would be violated by the decay of 
the J/tp to D*p. If the signal were due to J/ip decays with misidentified particles, the reconstructed mass would 
no longer lie at the J/tp value. 
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Figure 5: (a) Momentum distributions for all charged particles yielding M(D*p) values falling 
in the signal and side band regions of M(D*p), as defined in section 3.2, when combined 
with D* candidates of opposite charge, (b) M(D*p) distribution for p(p) > 2 GeV, with 
no proton dE/dx requirements. The data are compared with a two-component background 
model in which "wrong charge D" K ± tt ± combinations are used to describe non-charm related 
background and the ll D* MC" simulation describes background involving real D* mesons. 
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with no requirement on the proton likelihood. A strong signal is observed, with a reduced 
background which remains well described by the background model. The peak value and width 
of the observed signal are compatible with those for the standard selection shown in figure 2a. 

3.4 Photoproduction Analysis 

The analysis has also been carried out using an independent sample provided by HI data from 
the photoproduction region, where the scattered electron passes at a small angle into the back- 
ward beampipe, implying Q 2 < 1 GeV 2 . The hadronic final state is used to reconstruct y [22] 
and the selection 0.2 < y < 0.8 is imposed. The combinatorial background to the D* selec- 
tion is significantly larger for photoproduction than for DIS. To compensate for this, tighter 
proton and D* selections are imposed. The cut on the D* transverse momentum is modified to 
p T (D*) > 2 GeV. The region 1.6 < p(p) < 2.0 GeV, around the point at which the dE/dx 
parameterisations for protons and pions cross, is excluded and the requirement L p > 0.25 is 
made elsewhere. 

The distribution in M(D*p) for opposite-charge D*p combinations in photoproduction is 
shown in figure 6. Again, a clear signal is observed near M(D*p) = 3100 MeV, with mass 
and width compatible with those in the DIS case. The background distribution is reasonably 
modelled by the "wrong charge D" selection. The photoproduction signal is also separately 
observed forp(p) < 1.6 GeV and for p(p) > 2.0 GeV. 




Q I sH I I I I I I I I I I I I 

3 3.2 3.4 3.6 

M (D*p) [ GeV ] 

Figure 6: M(D*p) distribution from the photoproduction analysis, compared with a background 
model derived from "wrong charge D" K ± tt ± combinations. 
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4 Signal Assessment 



Fits to the M(D*p) distribution from the opposite-charge D*p combinations in DIS have been 
carried out to evaluate the peak position, width and statistical significance. Assuming the mea- 
sured width to be dominated by the experimental resolution, a Gaussian distribution is used 
for the signal, with the peak position, the width and the normalisation as free parameters. The 
background is parameterised with a power law of the form a [M(D*p] — m(D*)] 13 , with a and 
P as free parameters. A log-likelihood fit is made in the range 2950 < M(D*p) < 3600 MeV. 

The results of this fit are compared with the data in figure 7. They are also summarised 
in table 2, together with the results of separate fits to the D*~p and D* + p contributions. The 
fit yields a peak position of M(D*p) = 3099 ± 3 (stat.) MeV. The root-mean-square (RMS) 
width of the Gaussian is 12 ± 3 (stat.) MeV, compatible with the experimental resolution of 
7 ± 2 MeV, as determined from a simple simulation of the observed resonance with zero width 
and an isotropic decay distribution. The signal consists of N s = 50.6 ±11.2 events, from which 
the observed D*p resonance is estimated to contribute roughly 1% of the total D* production 
rate in the kinematic region studied. The fit results are not significantly affected when the 
background parameterisation is replaced with a polynomial or when the full distribution is fitted 
with the inclusion of a function to describe the rise at threshold. The results are stable against 
shifts in the binning, changes to the bin width in the range 1 MeV to 20 MeV and variations in 
the selection criteria which do not significantly alter the signal-to-background ratios for the D* 
or proton candidates. 



Sample 


Mass 


Gaussian width 


N s 




(MeV) 


(MeV) 




D*+p + D*~p 


3099 ± 3 


12 ±3 


50.6 ± 11.2 


D*~p 


3102 ± 3 


9±3 


25.8 ±7.1 


D*+p 


3096 ± 6 


13 ±6 


23.4 ±8.6 



Table 2: Results of the fit as described in the text for opposite-charge D*p combinations. The 
fitted position and Gaussian RMS width of the peak are given, together with the total number of 
signal events (N s ). The statistical uncertainties quoted take account of the correlations between 
the variable parameters in the fit. 

The systematic uncertainty on the mass of the peak is 5 MeV, estimated from the recon- 
structed masses of known states, such as the J/ip, with decays to particles in a similar momen- 
tum range, and from the variations in the peak position when the fitting procedure or selection 
criteria are modified. 

The probability that the background distribution fluctuates to produce the signal has been 
evaluated by comparing the observed number of events with background estimates for a win- 
dow spanning 3075 < M(D*p) < 3123 MeV, corresponding to ±2cr about the peak position 
according to the fit. The total number of events in this interval is 95. The background contribu- 
tion estimated using the fit described above is Nb = 45.0 ± 2.8 (stat.). A parameterisation of 
the background model shown in figure 2a yields a consistent value for N b . A more conserva- 
tive approach is to fit only the power-law background function to the full M(D*p) distribution. 
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3 3.2 3.4 3.6 

M(D*p) [ GeV ] 

Figure 7: M(D*p) distribution from opposite-charge D*p combinations in DIS, compared with 
the results of a fit in which both signal and background components are included (solid line) 
and with the results of a fit in which only the background component is included (dashed line). 

The data are compared with the result of such a fit in figure 7. The corresponding background 
estimate is N b = 51.7 ± 2.7 (stat.). The probability that a background of N b = 51.7 events 
fluctuates to produce at least the number of events in the signal is 4 • 1CT 8 , assuming Poisson 
statistics. This probability corresponds to 5.4 a when expressed as an equivalent number of 
Gaussian standard deviations. From the change in maximum log-likelihood A(ln C) when the 
full distribution is fitted under the null and signal hypotheses, corresponding to the two curves 
shown in figure 7, the statistical significance is estimated to be a/ —2 A(ln£) = 6.2 a. 

A state decaying strongly to D*~p must have baryon number +1 and charm —1 and thus 
has a minimal constituent quark composition of uuddc. The observed resonance is therefore a 
candidate for the charmed analogue 9® [5, 8] of the 9 + . The narrow width is reminiscent of that 
in the strange case. Given the relatively large mass of the resonance, it is also a candidate for 
an excited state such as the 6*° with spin 3/2 [23]. 
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5 Summary 



An investigation has been carried out of the invariant mass combinations of D* and proton 
candidates using HI deep inelastic electron-proton scattering data. A clear and narrow reso- 
nance is observed for both D* ~p and D* + p combinations with an invariant mass of M(D*p) = 
3099 ± 3 (stat.) ± 5 (syst.) MeV. The probability for the background distribution to fluctuate 
to produce a signal as large as that observed is less than 4 • 1CT 8 . The region of M(D*p) in 
which the signal is observed contains a richer yield of D* mesons and exhibits a harder proton 
candidate momentum distribution than is the case for side bands in M(D*p). The measured 
RMS width of the resonance is 12 ± 3 (stat.) MeV, consistent with the experimental resolution. 
A signal with compatible mass and width is also observed in an independent photoproduction 
data sample. 

The resonance is interpreted as an anti-charmed baryon decaying to D* ~p and its charge 
conjugate decaying to D* + p. The minimal constituent quark composition of such a baryon is 
uuddc, making it a candidate for a charmed pentaquark state. 
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